FACTORIZATION OF LAPLACE OPERATORS ON HIGHER 
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Abstract. This paper deals with the problem of factorizing integer powers 
of the Laplace operator acting on functions taking values in higher spin rep- 
resentations. This is a far-reaching generalization of the well-known fact that 
the square of the Dirac operator is equal to the Laplace operator. Using 
algebraic properties of projections of Stein- Weiss gradients, i.e. generalized 
Rarita-Schwinger and twistor operators, we give a sharp upper bound on the 
order of polyharmonicity for functions with values in a given representation 
with half-integral highest weight. 



1. Introduction 

Clifford analysis provides a generalization of complex analysis in the plane to a 
higher-dimensional setting, in which the role of the Cauchy-Riemann operator is 
played by the Dirac operator. It is centered around the study of functions on the 
vector space M™ taking values in its Clifford algebra or the corresponding spinor 
representation, see e.g. [Tl [71 for the standard references. Let {ei, . . . , Cm} denote 
an orthonormal basis for the Euclidean vector space M™. The Clifford algebra 
Cm of the complexified vector space C™ is generated by the multiplicative relations 
CiCj + CjCi = —2Sij and has the structure of a graded vector space C„i = ®5!Lo ' 

(k) 

where Cm is spanned by elements e^^ . . . e^j. with 1 < ii < . . . < ik < m. The space 

(2) 

Cm , provided with the commutator bracket, defines a model for the Lie algebra 
so{m) of the group Spin(r7i), which can be defined as the Lie group consisting of 
all even Clifford products of unit vectors in M™. The spinor representation S of 
Spin(m) can be realized as a minimal left ideal in Cm by means of a standard con- 
struction involving the Witt basis. For m odd, § is irreducible with highest weight 
. . . , i), for m even it decomposes into two irreducible summands S+ © with 
weights (^, . . . ± . For the sake of convenience, we will restrict to the case of odd 
dimension in this paper (although it should be pointed out that the results also 
apply to the case of even dimensions, taking parity changes into account). 
The Dirac operator :— J27Li^idxi acts on the space C°°(R™,§) of §- valued 
functions as a Spin(TO)-invariant elliptic differential operator, and its kernel is de- 
fined as the space of monogenic functions. In view of the fact that = —Am, the 
Dirac operator factorizes the Laplace operator in m dimensions, which means that 
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monogenic functions are also harmonic. 

In recent years, several authors O IH El [Til El [9l [8] have been working on general- 
izations of Clifford analysis techniques to the so-called higher spin representations: 
these are arbitrary irreducible representations S;^ of Spin(m) with dominant half- 
integral highest weight 

where A contains integers only. Given this highest weight, there is (up to normal- 
ization) a unique invariant first-order differential operator 

■Rx :C°°(R™,Sa) ^C°°(M'",§a) • 

It is given as the projection of the Stein- Weiss gradient of a function to the irre- 
ducible summand Sa in (M'")*«)Sa. For A = (1, 0, . . . , 0) := (1) the operator 11(^1) is 
the multidimensional analogue of the classical Rarita-Schwinger operator, see [RS]. 
In full generality, the operator TZ\ is called the higher spin Dirac (HSD) operator 
for the representation §a of Spin(m) or its Lie algebra i3„, see Section 3. 
The main result of this paper is a formula for the factorization of a suitable power 
of the Laplace operator on an arbitrary higher spin representation through HSD 
operators: 

Theorem 1. Let X be a dominant integral weight, and let n £ N be an integer 
satisfying n > Xi. Then there exists an invariant differential operator A\^n acting 
on C°°(M",§a) such that 

A" = 7^AAA,n7^A . 

Thus every function in KeiTZx is (Ai + l)-polyharmonic. The structure of KerT^A 
is exphcitly known only for the cases A = (k), see e.g. [4], and A = {k,l), see e.g. 

In both cases the bound (Ai -I- 1) is optimal, i.e. there are functions in KerT^A 
which are not Ai-polyharmonic. 

Our paper is structured as follows: we will briefly explain how the HSD operators 
arise within Clifford analysis framework in Section 2. At the same time, we will 
show how knowledge on the factorization of the Laplace operator on the underlying 
higher spin representations can be used to describe the (polynomial) homogeneous 
kernel spaces. Sections 3 and 4 are then devoted to the actual factorization of 
the Laplace operator. To do this, we will make use of techniques coming from 
representation theory. Therefore, we will briefly fix the notations in Section 3, 
before turning to the main problem in Section 4. 

2. Higher spin Clifford analysis 

Let <Cm be the universal complex Clifford algebra generated by an orthonormal 
basis {cp : 1 < p < m} for the vector space R"*, satisfying the multiplicative rela- 
tions CpBg -|- CqCp — —26pq. As was pointed out earlier, the spinor space § can then 
be realized as a minimal left ideal S — Cm/, where / denotes a primitive idempotent 
(seeH). 

The main motivation for considering HSD operators within the framework of Clif- 
ford analysis comes from the crucial fact that any higher spin representation §a 
can be explicitly realized as a vector space of §-valued polynomials in several vec- 
tor variables Up G M™, satisfying certain systems of equations. In what follows, the 
variables Up and their associated (classical) Dirac operators dp := 9„ have to be 
seen as dummy variables, "fixing the values" . 
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Definition 1. A function f : M*^™ — > S, {ui, . . . ,Uf.) ^ /("i, • • • , Wfe) is simplicial 
monogenic iff the following conditions are satisfied: 



The vector space of S-valued simplicial monogenic polynomials which are homo- 
geneous of order Xp in the variable will he denoted by Sa — ^Xi_,...,Xk (with 
Ai > . . . > Afc > from now on). 

In [CSV], the following result was obtained: 

Theorem 2. With respect to the regular representation L of the spin group on Sx, 
given by L{s)f(ui, . . . , Uj.) — sf{sUiS, . . . , sw^s), the space Sa defines a model for 
the the irreducible Spin{m) -representation with highest weight A'. 

This means that a HSD operator in Clifford analysis can be seen as an operator 
TZ\ acting on functions /(ai; Wi, ■ • ■ iMfc) £ S\, such that 



with p and q indices as in the definition. Let us then first consider the (special) 
cases mentioned earlier, i.e. for which the kernel spaces are fairly understood. 

2.1. The Rarita-Schwinger operators TZk- Fixing a highest weight A' = (k)' , it 
follows from the Theorem above that the irreducible module with highest weight A' 
can be realized as the space Mk 'Pfe(K™, S) nker 9„ of fc-homogeneous monogen- 
ics. The corresponding HSD operator, generalizing the Rarita-Schwinger operator 
coming from physics |10| . is usually denoted as TZk and is then defined as 



The function theory for this operator was to a great extent developed in [3j |4] , in 
which the fundamental solution, the Cauchy integral formula and explicit descrip- 
tions for polynomial solutions were obtained. In particular the following result was 
obtained: 

Theorem 3. There exists a differential operator ^2fe+i acting between spaces of 
A4k-vo,lusd functions, such that TZkA2k+i = ^^n^ ■ 

In order to illustrate the link between this result and the structure of the space 
kerji TZk of ft,- homogeneous (polynomial) solutions for the operator TZk , we also 
mention the following crucial result: 

Theorem 4. The induction principle states that 



kei'hTZk ^ Mh,k ® d^^iui^eih^iTZk^i) ■ 
Here, the Aih,k stands for the space of {h,k) -homogeneous polynomials which are 



dpf 

{Up,dq)f 








l<p<k 
i < P < q < k 



dp{TZ\f{x;Up)) = = {Up,dg){TZxf{x;Up)) , 





Remarks: 
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(1) This theorem exhibits an important feature of HSD operators, as opposed 
to the classical Dirac operator: the space of polynomial solutions is no 
longer irreducible under the action of the spin group. As a matter of fact, 
decomposing this kernel space is a central problem within this setting. 

(2) The action of the operator d^^ from the Theorem is nothing but the action 
of a (dual) twistor operator, cfr. infra. 

(3) In particular, the theorem above illustrates that solutions for TZk are indeed 
polyharmonic: it suffices to perform induction on the index k to see that 
A^+-'^ indeed acts trivially on kerTZk- 

The Theorem essentially says that in order to decompose the kernel of the oper- 
ator TZk, one needs two basic ingredients: a special type of solutions (the double 
monogenics), together with information on the action of the twistor operator on 
solutions for the "simpler" system TZk^if = 0. 

2.2. The higher spin operators TZk.i- In a recent series of papers [5J[n], we have 
considered the case where A' — {k, I)' . Our motivation for doing so lies in the fact 
that the case of the Rarita-Schwinger operator, does not yet contain the germs from 
which the behaviour of the most general case can be conjectured. The corresponding 
HSD operator TZk^i, acting on the space of iSfc^^-valued functions taking, was found 
to be given by: 

The main function theoretical results were developed, and in particular an accurate 
description of the polynomial solution space was obtained. Without going into full 
details, we mention the following result: 

Theorem 5. A refined version of the induction principle leads to the following 
decomposition for the kernel of the operator TZk.i '■ 

ker,7^fe,^ ^00^^^,_^.,_, . 

1=0 j=0 

For {i,j) = (0,0), one obtains the space A^^j,/, which is the analogue of the 
space of double monogenics in the present setting. To be more precise, this is 
the space of triple monogenics, satisfying an additional constraint (hence the extra 
superscript) to make sure that the correct values are assumed. For all other choices 
{hi) 7^ (0, 0), one obtains solutions for "simpler" equations TZk-i,i-jf = which are 
in certain sense inverted through a suitable action of associated twistor operators. 
It is then crucial to point out that not all the indices (k — i,l — j) are used in the 
Theorem above. As a matter of fact, the indices needed in the decomposition above 
are precisely the indices needed in the decomposition of the Laplace operator on 
5fe,;-valued functions (see Section 4 for details). 

2.3. The general case TZ\. In a recent paper [6j, techniques were developed to 
obtain explicit expressions for general HSD operators Ti\. Moreover, the authors 
obtained the corresponding special class of solutions M'',^ ^ for these operators which 
generalizes the space of double monogenics from the Rarita-Schwinger case. To 
complete the description of the polynomial kernel space, one then only needs to 
understand the behaviour of the twistor inversion, relating solutions for TZ\ to 
solutions for "simpler" HSD operators. A crucial step towards this understanding 
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is the Theorem we are about to prove in Section 4, regarding the decomposition of 
the Laplace operator. 

3. Notations 

Let be a complex semisimple Lie algebra of the type i3„ with Cartan subalgebra 
[). By (e^)" wc denote the standard basis of C" = f)*. The dominant part A^; of 
the weight lattice A of g can be expressed in the standard basis as J2i ^i^i '^ith 
Xi either all integral or all half-integral non-negative numbers satisfying Ai > A^+i 
for alH e {1, . . . ,n — 1}. For A integral, this is the highest weight of a tensorial 
representation with symmetry given by the Young diagram with A, boxes in the 
i-th row. Let us denote by A' the weight A + i e. the Cartan product of A 

and the spinor representation. Finally, we denote by Va (resp. S^) the irreducible 
representation with highest weight A (resp. A') and by V\ (resp. Sx) the set of 
smooth functions with values in this representation. 

Next, we also define a graph W whose set of nodes consists of all integral dom- 
inant weights and where two weights A,/i are joined by an edge iff '^"^li^i — 
= 1- The graph VV has the same set of nodes, but the edges are defined by 
S^r=i('^» ~ '"»)^ < 1. It is then easily seen that there exists a first-order 0-invariant 
operator from Va to iff A and /i arc joined in W, and wc denote such operators 
by -Da.m : — >■ Va- Each of these operators is unique, up to a constant: the issue 
of normalization in the cases of interest will be dealt with later. As we are mainly 
interested in operators between half-integral representations, generalizing twistor 
and Rarita-Schwinger operators, wc introduce two more symbols: 

Tx.^i ■■= Dy,f,' and Tlx := Dx',x' ■ 

We also introduce a similar symbol A a for the Laplace operator on Sx, and for the 
sake of convenience, we put any operator Tx,ij,,Rx,^\ for which the weights are 
non-dominant equal to zero. 

4. Path operators 

First of all, we define two particular sets of weights: 

Definition 2. Let ii and v he two integral dominant weights such that fj, y u in 
the Bruhat order, i.e. j-ii > Vi for all i G A path from v to is a 

sequence {Xp)q of weights in W such that Aq = /i. Xn = v, and for all indices 
p € {1,...,A'^}, we have that Xp >- Xp-i and jAp.Ap_i| = 1. This means that 
for two successive weights Ap_i and Xp, there is an index Ch{p) S {1, . . . , n} such 
that Xp — Xp-i — tch(p), we call this index the change at p. A reverse path is a 
sequence of weights {Xp)^ such that (Ap)^ is a path. When there will be no danger 
of confusion, we will use the term, path for both paths and reverse paths. We call 
the number of edges N the length of the path or reverse path and denote it by 
clearly it is the distance between fi and v in the Manhattan metric for any path 
between fj, and v. 

Definition 3. Given a dominant integral weight /i. vje define the box B(fi) of the 
weight fi as the set of dominant integral weights X satisfying Hi > Xi > /Ltj+i for all 
i € {1, . . . ,n — 1} and > Xn > 0. 

Given a path, we then define the following differential operator: 
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Definition 4. For any path {Xp)o , the path operator is a composition of twistor 
operators along the path: 



Similarly, the operator Gi^^^ is defined by composing along the reverse path. 

The notation from the definition above is justified by the following crucial 
lemma: 

Lemma 1. Let ^ > v he dominant integral weights. Then there exists a normal- 
ization of the twistor operators Tx^ui such that 

(i) and G^^fj, do not depend on the choice of the path from v to \x, modulo 

a constant multiple. 

(ii) For a suitable choice of normalization of the twistor operators, one has that 



(ill) For all integral dominant weights X -< fi not in B{^), the path operators 
Gij.,x and Gx^^ are zero. 

Remark: note that property {ii) from above was already mentioned in a previous 
section. It essentially relations solutions for different HSD operators, through the 
action of a twistor operator. 

Proof. The tool for proving these properties is the decomposition of the twisted 
Dirac operator. For any dominant (integral) weight A, the representation Yx®^ 
decomposes as 0V^, the sum being taken over the set 



where each summand appears exactly once, up to possible omissions of weights 
which are not dominant. For a given dominant weight A, let us encode summands in 
this decomposition into a sequence of pluses and minuses like (— , +, +, — , — , — , +). 
We then consider the operator D = id® onNx® §-valued functions, where 
is the Dirac operator. This twisted Dirac operator D is then equal to a sum of 
twistor operators (and HSD operators) on each summand. These are nonzero iff the 
plus-minus codes of the source and target summands differ at most at one position. 
If we think of these 2" (±)-codes as vertices of an n-dimcnsional hypercube, there 
is a nonzero twistor operator between two different vertices iff they are joined by an 
edge of the hypercube. Since = —A is a scalar operator which preserves each 
summand, we have all crucial information encoded in diagrams like this: 




iiV-2 



'0 





(+,+) 



(-,+) 




?(+,-) 
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Each square with nodes (++),(H — ), ( — h) and ( ) hereby represents the n- 

dimensional hypercube of summands and the arrows describe aU nonzero twistor 
(and HSD) operators in the decomposition of the twisted Dirac operator. 
For K, t G Aa let denote the restriction to the summand G Va CE> § and p„ the 
projection to in the same diagram. If we then choose the normahzation of the 
operator T„_t in such a way that r„ = o D o r^, the scalar operator p„ o o 
gives rise to three operator identities. First of all, for k = t we get that 

(1) - A, = 7^2 + T^:^T^.,, . 

K,Cj| — 1 

There are precisely n terms in the sum, but some of them can be zero (when uj is 
not dominant). Secondly, for \h,l\ = 1 we get 

(2) T^,ji, + 7^«^«,, = . 

Finally, for |k, = 2 we have that 

(3) T^,eTe,L + TK,.t^T^.L = , 

with 6 and oj the only two weights in on a distance 1 from both k and l. 
Summands with |k, i| > 2 do not give further relations. 

Let us consider a specific path (A^)^ between /i and >~ going through the 
weights: 

Aq = {vi,V2,V'i,...,Un-l,yn) 
Aji = (/il, 1^2, ^^3, ■ • ■ : Vn-l, Vn) 
Ai2 — (/^l; /^2, V^, ■ ■ ■ , Vn-l, Vn) 



A'i„_i — (MIi M2, Ma, ■ ■ ■ ^ fJ-n-1: ^n) 
Aat = (/^l, ^2, M3, ■ • • , 

where = io < ii < *2 < • ■ - in-i ^ in = N . These weights determine the path 
uniquely, since between any successive pair of them there is a unique path. The 
sequence {Ch{i)}^ is non-decreasing and increases at each ik, k € {1, . . .,n — 1}. 
For a path (i^i)^ from i/ to n different from the given path (A,)^, there is always 
an index j £ {1, . . . , N} such that Ch(j) < Ch{j — 1). The path (cji)^ can then 
be deformed, by changing 

ujj^i = ujj^2 + ech{j-i) = - ec/iO) 
into a new weight 9, given by 

6 := Wj_2 + ech{j) = ^3 - ec/iO-i) ■ 
Using operator identity ([3]), we then have that 

-L UlM ,i^N -1-'- i^N -l^t^N -2 ■ ■ ■ ^ t^j ,^^3-1-'- i^j-l.^j-l ■ ■ ■ ^UJl;U]o 

equals 

In such a way it is clearly possible to deform any path from v to fji into (Ai)^ and 
the resulting path operators differ only by a multiple coming from the signs and 
the way how different twistor operators are normalized. An analogous argument 
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works for path operators defined by a reverse path. 

We will now choose the normalization of the twistor operators in such a way that 
G^^u will no longer depend on the path. Let us therefore write a weight /i as 
{^i, ^2, ■ ■ ■ , ^J'm, 0, . . . , 0), with > for m G {!,...,«}. We then define G^,o as 
the composition of twistor operators along the unique path through 

Ao = (0,0,0,. ..,0) 

A,, = (/ii,0,0,...,0) 

Ai2 = (mi,M2,0,...,0) 

Aat = {Hl, • . ■ , fJ-m, 0, . . . , 0) , 

and choose an arbitrary normalization for G^.o along this preferred path. For any 
index p < m, we can now normalize the twistor operator T^+e^.^ by the requirement 
Gfi+ep.o — Tji+epiM^M.o- This makes the complete diagram of twistor operators 
commutative and also chooses the normalization for any path operator G^^^. The 
same way we can normalize path operators for reverse paths. 
If (A)^ is a path from i' to fi then by the relation ([2]) we get that 

and continuing this process we establish the second assertion up to a multiple. An 
arbitrary normalization of TZq then gives us a consistent choice of overall normal- 
ization of other operators TZ^ by requiring that 

Ti-fiG^fl = G^flTZo . 

To prove the third assertion, we need to show that for any A ^ /i such that A ^ B{ii), 
the operator G^^a can be expressed via a path in which a composition of two specific 
operators is zero. If A ^ B{fi) then for some i £ {1, . . . , n—1} one has that A^ < /i^+i 
and Hi+i > fJ.i+2 (with the additional definition jin+i 0). We choose the path 
from A to ^ as a succession of paths from A to A_, then the path (A_, Aq, A-|_) and 
then from A+ to fi, where 

A- = (/il, fi2, ■ ■ ■ , Mi-l, — 1, fJ-i+l — 1, • • • , fJ-n-l, fJ'n) 

Aq = (M1, M2, ■ • ■ , ~ 1; ■ • ■ , Mn-1, /^n) 

A+ = (/il, /i2, ■ • • , fJ-i+l, IJ-i+2, • • • , fJ-71-1, l^n)- 

We see that all three weights are dominant and Aq G B{ii), since > /i^+i > Mt+i, 
Mi+i ^ Mi+i ~ 1 ^ and the other coordinates are the same as in fj, itself. 

Similarly A+ € B{^). As A_ :^ A and fi >- A+, the path operators Ga_,a, G^,a+ are 
well defined and 

Gp,A = G^^a+Ta+.Ao^^Aq.A- Ga_,a 

Now a weight 

K (/il - l,/i2 - l,...,/ii-l - - 1,/it+l - 1, . . . ,/i„-l,/i„) 

is still dominant and the weights A_, Aq, A+ are sums of k with 

= (1,1, ...,1,0,0,0,. ..,0,0) 
^0 = (1,1, •••,1,1,0,0,. ..,0,0) 
UJ+ = (1,1, ...,1,1, 1,0,. ..,0,0), 
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respectively, i.e. the representations §a_ , Saq i Sa+ are Cartan products of with 
the representations 
operators 



cj_ , VcjQ I ■ The sequence of first order invariant differential 



is in fact part of the de Rham complex and so the composition is zero. If we twist 
this sequence by S^j the composition operator is still zero and so is its restriction 
and projection to the Cartan parts S\_ and respectively. Since the only weight 
^ Ad such that 16*, A+| — 1 and |^?, A_| = 1 is Aq, the twisted de Rham operator 
restricted on S\_ an projected to Sx^ must be a composition of twistor operators 



Sx 



Sx 



which is up to a multiple Tx^,AoTao,a_ • Hence G^,a = 0. A similar argument shows 



that Gx.^, ^ 0. 

Theorem 6. Let /i G A^, then for n> fii 



□ 



A" 



J2 c(m,A)G^,aA; 

AeB(At) 



I— I Ai, A| — 1 



G 



7^. 



where c{pL, A) are constants. 

Proof. Relation ([1]) for the twisted Dirac operator on ® § gives: 

X-<fj, 
Ai,A| = l 

where we have to include constants c{fj,, n) and ci{fi,X) to compensate for our 
choice of normalization of the RS and twistor operators, and for the sign included 
in relation ([Ij . If we expand the Laplace operators in the terms of the summation 
over A in a similar way, we get: 

+ J2 c{ii,\)G,,^xnxAl-^nxGx.^.+ C2{ii,\)G^,^xAl-^Gx.^. 



X-<ti 
|ai,A| = 1 



A^A^ 
Ai,A|=2 



We can now use the second part of Lemma [T] to switch the HSD and path operators 
in the first sum and we can then expand the Laplace operators in the second sum 
again. If we continue plugging in and expanding, we finally arrive at 



A" — 7? 



E 

X^fi 
Ai:A|<r! 



c(A*,A)G^,AAr'''^^'"'GA.A. 



E 

X-ifj. 
|Ai,A|— n 



c„(/i, A)G^,aGa„ 



The second sum is zero since the element (fi2, /^3, • ■ • , Mm 0), which lies at maximal 
distance from fi G B{p), has distance and n > fii. A similar argument allows 
us to keep only the terms corresponding to A G B{^) in the first sum. □ 

As a result, we have the following conclusion, which generalizes the fact that all 
monogenic functions in the kernel of are harmonic: 
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Corollary 1. Given a dominant integral weight /i G and a function f{x) £ Sx, 
we have that 

nxf{x) = ^ A^i+i/(x) = . 

Let us recall that for the cases of A = (A:) and A = {k,l), where the kernel of 
TZ\ is known explicitly, such polyharmonicity result is a direct consequence of the 
decomposition of the kernel into a direct sum of special solutions for "simpler" HSD 
operators TZf^. The set of highest weights /x that take part in this decomposition is 
precisely B{\), i.e. it is the same set that appears in the corresponding factorization 
of the operator A^. This suggests how the structure of the kernel of TZx may look 
like in general, and indeed, dimension checking in randomly chosen cases supports 
this conjecture. We will treat this in future work. 
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